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Abstract 



I/-} , We show that the strong coupling limit of d-dimensional quantum electrodynamics 



with 2 d /2t d / 2 l flavors of fermions can be mapped onto the s=l/2 quantum Heisenberg 
antiferromagnet in d-1 space dimensions. The staggered Neel order parameter is the ex- 



the flux phase and the conventional Neel ordered phase of an antiferromagnetic t-J model. 
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pectation value of a mass operator in QED and the spin-waves are pions. We speculate 
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that the chiral symmetry breaking phase transition corresponds to a transition between 
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The possibility that quantum electrodynamics (QED) can have a nonperturbative 
ultraviolet fixed point has been investigated by many authors [1-8]. Such a fixed point 
would give QED sensible ultraviolet behavior and save it from triviality by avoiding the 
Landau ghost (or Moscow zero) [9-10]. The most popular candidate for an ultraviolet fixed 
point arises from a second order phase transition which occurs when the electromagnetic 
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coupling is increased to a critical value where the theory breaks chiral symmetry. 

The physical mechanism for such a phase transition advocated by Miranksy [3] is 
that of 'collapse' of the electron-positron wave-function when the charge reaches a super- 
critical value, e 2 ~ 4w. There are two ways to screen super-critical charges. In the case of 
a super-critical nucleus the high electric field produces electron-positron pairs, ejects the 
positron and absorbs the electron to screen its charge. Alternatively, the pair production 
is suppressed by fermion mass, so the system can stabilize itself by increasing the electron 
mass, thus the tendency to break chiral symmetry. These ideas are supported by studies 
of the Schwinger-Dyson equations for QED in the quenched ladder approximation [3-5] 
which find a line of fixed points in the e 2 -plane between and e 2 = 4-7r 2 /3 where the 
theory breaks chiral symmetry dynamically and has a nontrivial continuum limit. This 
critical point has other interesting behavior such as large negative anomalous dimensions 
for fermion composite operators so that some four- fermion operators are relevant [4] . There 
is also some numerical evidence for this phase transition [7]. 

In this Letter we shall present a lattice model whose weak coupling continuum limit is 
4-flavor QED with light fermions and which, in its strong coupling limit, exhibits sponta- 
neously broken chiral symmetry. Conventional analysis of Euclidean lattice gauge theory 
formulates it as a classical statistical mechanics problem and seeks second order phase tran- 
sitions so that a nontrivial continuum limit exists [11]. Here we formulate a Hamiltonian 
version of lattice QED as quantum statistical mechanics [12]. In the latter a relativistic 
continuum limit exists if there are gapless degrees of freedom and if those degrees of free- 
dom have a relativistic dispersion relation, uj(k) ~ \k\. This is, of course, true for the weak 
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coupling continuum limit of the lattice model which produces QED with 4 flavors of elec- 
tron (reminiscent of the standard model with four generations). We show that the strong 
coupling limit is equivalent to the spin 1/2 quantum Heisenberg antiferromagnet and that 
the Neel order coincides with chiral symmetry breaking. In three or greater dimensions, 
there is a rigorous proof that the ground state of the s=l/2 antiferromagnet has Neel 
order [13]. Also, there is evidence from numerical simulations and large spin expansions 
for Neel order in two dimensions [14] . Furthermore, these systems exist in nature and can 
be studied by experiment. In dimension 3 and higher the order persists for some range of 
temperature. 

This result indicates that in strong coupling 4-flavor QED has a nontrivial continuum 
limit. The light excitations are the 'pions' which coincide with the spin waves of the 
antiferromagnet. They have a relativistic dispersion relation and interactions which are 
commonly represented by a nonlinear sigma model. Corrections to the strong coupling 
limit take into account fermion hopping terms similar to those in a gauge invariant t-J 
model. We conjecture that if 1/e 2 is increased to some critical value the chiral symmetry 
of QED is restored. The resulting phase of the t-J model is known as the flux phase which 
is a plaquette-centered antiferromagnetic gapless semiconductor rather than site-centered 
antiferromagnetic insulator. 

We shall use staggered fermions on a (d-l)-dimensional lattice and continuum time 
which are obtained by spin-diagonalization [15] of the naively latticized Dirac Hamiltonian 

Hf = — (^\x)ot? "Vjtfj(x) — V J (x)^ 'ip(x)) = — — (ip\x + j)a^ip(x) — ifi \x)a? ifi(x + j) 

x,j x,j 

(1) 
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where V J is the forward lattice difference operator, Vjf(x) = f(x + j) — f(x), a J are the 
2 [d/2] x 2l d / 2 l Hermitean Dirac matrices, x, y, . . . refer to sites on a hypercubic lattice, and 
. . . refer to unit vectors. (Here [d/2] is the largest integer < d/2.) The second form 
of the Hamiltonian in (1) describes a fermion hopping problem in a U(2t d / 2 ]) background 
gauge field given by the unitary matrices a J . The crucial observation which allows spin 
diagonalization is that this background field has only U(l) curvature, i.e. if we consider 
the product around any plaquette a J a k {a^)^ (a k Y = —1. This allows diagonalization 
using the gauge transformation ip(x) — > (a 1 ) 11 (a 2 ) X2 . . . (a < ^ d ~ 1 ^) X(d - 1 '>i(j(x) resulting in the 
Hamiltonian 

H f = ~\ E(" 1 ) Erl Xp (V 0* + 3>W - ^(^(x+j)) (2) 

which describes 2^1^ species of lattice fermions with background U(l) magnetic flux n 
through every plaquette of the lattice. Each species of fermion must have the same spec- 
trum as the original one given by the Dirac Hamiltonian (1). This allows reduction of the 
fermion multiplicity by a factor of 2^ d / 2 \ The result resembles a condensed matter hopping 
problem with a single species of fermion where there is a background magnetic field 7r per 
plaquette. 

Chiral symmetries are obtained by lattice translations by one site. This translation 
interchanges the even x p =eveii) and odd x y =odd) sublattices. The substi- 

tutions 

ip(x) -> (-1)^+1 Xp ^{x + j) (3) 

leaves the Hamiltonian in (2) invariant. A candidate for Dirac mass operator, which 
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changes sign under the transformations in (3), is the staggered charge density operator 

^ = E(- 1 ) E " lx ^ t W(^) (4) 

X 

If we introduce N species of lattice fermions, the continuum limit of (2) describes 2 d ~ 1 N/2^ d ^ 2 ^ 
species of massless Dirac fermions. 

To obtain the continuum limit and the number of fermion species, we first divide the 
lattice into 2 d ~ 1 sublattices according to whether the components of their coordinates are 
even or odd. For example, when (d-l)=3, we label 8 fermion species as ip(even, even, even) = 
fa, •(/'(even, odd, odd) = fa, '(/'(odd, even, odd) = fa, ^(odd, odd, even) = fa, -0(even, even, odd) 
-05, if; (even, odd, even) = fa, if; (odd, even, even) = fa, ip (odd, odd, odd) = fa. Then, if we 
add the mass operator in (4), in momentum space the Hamiltonian has the form 

H f = [ d 3 k if)\k) (A l sink l + Bm)ifj(k) (5) 

a*\ (1 

a 1 )' I -1 

a 3 

a 1 are Pauli matrices, we have used the Fourier transform fax) = J Qb pt-y*/ 2 e ~ lkx ^(k) 
and Qb is the Brillouin zone of the (even,even,even) sublattice, — n/2 < h < n/2. The 



a 1 \ ( 1 

where the 8x8 Dirac matrices are ip(k) = (fa, . . . , fa), A 1 = I t \ , B = 



= nl . = V\ , « 3 =( a « Is) («) 
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fermion spectrum is u(k) = sm2 ^ + m2 an d only the region ki ~ is relevant to 

the continuum limit. We have normalized if;(k) so that 

{fax),^(y)} = S(x-y) , {fak),^(l)} = 8(k-l) (7) 

If we define (3 = ^ ^ ^ 2 ^ an d the unitary matrix M = | ^ j ^ i |# ^ anc ^ ^ = 
Mfa with fa = [ip a , fa) the Hamiltonian is 



In the low momentum limit, sin&j ~ ki, with fermion density 1/2 per site, we obtain 2 
continuum Dirac fermions. Furthermore, the staggered charge operator gives a Dirac mass 
of differing sign for the two species. In d dimensions the proceedure is similar to this. In 
general, we shall consider N lattice species in d dimensions which yeilds 2( d_1 ) N/2^- d / 2 ^ 
continuum species of Dirac fermions where the lattice fermion density is N/2 per site. 

In lattice electrodynamics the gauge field Ai(x) and electric field E^x) associated 
with the link between x and x + i are conjugate variables, [^4j(x), Ei(y)\ = ie 2 5ij5(x — y). 
We shall represent this commutator by taking the quantum states as functions of Ai(x) 
and Ei(x) = —ie 2 d/dAi(x). The anticommutator algebra (7) is represented by a 2 iV -level 
fermion system at each site with cyclic vector defined by ip(x)\0 >= Vx. The Hamiltonian 
is 



H = Y.-\dX7xf + £ 2?*3 + £ + JY AAX) ^) + h-c) (9) 

(V. 1 3 \ / rp n A rp A 



where = ViAj—VjAi for noncompact QED, F^ = sm(V iA 3 ■— VjAA for compact QED, 
and t Xj i contains a background field of it (mod 2%) through every plaquette of the lattice. 
The Hamiltonian is invariant under the gauge transformation Ai{x) — > Ai(x) + Vix(x), 



Gauge invariance is imposed as a physical state condition, Q(x)^ p h ys (A) = 0. 

We shall first consider the case N = 2. Then there are 4 flavors of 4 component 
fermions in 4 dimensions and 2 flavors of 4 component fermions in 2+1 dimensions. For 
strong coupling perturbation theory we write the Hamiltonian as, H = Hq + Hi + Hi 



ip{x) — > e tx ^t(j(x) which is generated by 




(10) 



where I = -yE #i = EWH + h.c. and H 2 = ^ E^g 



are each gauge invariant 
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operators. The leading order ground state is the ground state of H which is A-independent 
and is gauge invariant when each state has fermion occupation number one, l^oba;]) = 
Y\ x ip Tx (x)\0 >. This is a normalizable state in compact QED where Ai(x) is integrated 
from to 2tt and is a non-normalizable component of a continuum spectrum in non-compact 
QED. It is highly degenerate: each of the 2 V (where V is volume) components labelled 
by [t x ] has the same eigenvalue of H . The degeneracy is resolved by diagonalizing the 
matrix elements of perturbations in 1/e 2 . First order perturbations to the vacuum energy 
vanish. Thus, the leading term in the vacuum energy is of order 1/e 2 , and is given by the 
lowest eigenvalue of the matrix 

5 2 E = -(Mrx}\H! 1 ffilfroK]) + (Mr,]\H 2 \Mr x ]) (11) 

-no — 

The second term is diagonal and therefore is irrelevent to resolving the degeneracy. Diag- 
onalizing the matrix in the first term is equivalent to solving the eigenvalue problem for 
the four-fermion Hamiltonian 

HeS = ~^Y1 I*<mI V(z + i)^{x)^{x)^{x + i) (12) 

x,i 

restricted to the subspace of the Hilbert space where each site is singly occupied. With 
the identity — 2ip^ (x)ip(y)ip^ (y)ijj(x) = ip^ (x)atp(x) ■ ijj^ (y)aip(y) + tp^ (x)il>(x)^ (y)ip(y) we 
obtain the Hamiltonian of the s=l/2 quantum antiferromagnet, 

2 . 

^eff = — ^2 l^,i| 2 V ,t ( a; + + i) ■ ^{x)aip{x) + const. (13) 

x,i 

Fermion mass operators are staggered charge density operators. For N=2, consider 
a mass operator with differing signs for the two lattice flavors, i.e. m J2 X (— l)^ x ^a s i(j 
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which in the naive continuum limit corresponds to the mass matrix 



— m 

V 



(Note that there is no chiral anomaly in this channel.) This corresponds to a staggered 
magnetization operator and its expectation value is the order parameter of the antiferro- 
magnet. There is a proof that for s=l/2 and (d-l)> 3, (and it is widely believed to also 
hold in d-l=2 although there is no rigorous proof) that the Heisenberg antiferromagnet 
has a Neel ordered ground state, i.e. 



In quantum electrodynamics, this is the chiral limit and implies that there is a nontrivial 
expectation value of the mass operator and thus spontaneous chiral symmetry breaking. 
In (d-l)> 3 this persists for some finite range of temperature whereas in (d-l)=2 it can 
only be true at zero temperature. 



breaks chiral symmetry. The only light particles in the spectrum are spin- waves which are 
the Goldsone bosons for breaking of the SU(2) symmetry of the Heisenberg model. In QED 
they are the 'pions' corresponding to the breaking of the chiral symmetry and electrons 
and all other charged excitations are confined. There are not enough spin waves to account 
for the breaking of SUl{4) x SUr(4) chiral symmetry. The reason for this mismatch of the 
number of Goldstone bosons and broken symmetries arises from the presence of relevant 
operators which reduce some of the apparent continuous chiral symmetries of the naive 
continuum theory to discrete ones of the lattice. These discrete symmetries are broken 
but do not require Goldstone bosons. The only true flavor symmetry of the lattice theory 



lim lim < 




(14) 



X 



We conclude that, in the infinite coupling limit, 2 d /2^- d ^ 2 ^ flavor (4 flavor in d=4) QED 
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is SU(2). The apparent SUl{4) x SUr{A) symmetry of the weak coupling limit stems 
from the fact that all of the operators which could break SUl(A) x SUr(A) to SU(2) plus 
discrete chiral symmetries are irrelevant in that limit. However, they can become relevant 
at strong coupling and the full SUl(A) x SUr(A) symmetry is absent. Also, recall that 
Goldstone's theorem applies to the continuous symmetries of the bare Hamiltonian. The 
effective Hamiltonian (13) is invariant under the discrete chiral transformations in (3). 
However, the ordered ground state is not - therefore the discrete chiral symmetries which 
involve translations by one lattice site are also broken as 1/e 2 — > 0. 

We can think of the higher order in 1/e 2 corrections as perturbations of the Heisen- 
berg model. When we increase 1/e 2 we eventually should arrive at a critical coupling 
where the chiral symmetry is restored. A future technical problem will be to estimate the 
critical coupling. Also, spin-wave analysis which is exceptionally good for the quantum 
antiferromagnet can also be used to look at the dynamics of strong-coupling QED. 

We expect that chiral symmetry breaking should persist for N in the vicinity of N = 2. 
However, when N=l (and whenever N is odd), since operator ip^ip has integer eigenvalues, 
it is not possible to find states which are annihilated by the gauge generator in (10) without 
nontrivial electric fields. The strong coupling ground state energy is necessarily of order 
e 2 and the ground state is the lowest eigenstate of the effective coulomb Hamiltonian, 

H ° = E y - 1/2)^^ - 1/2) (15) 

X 

We expect that the ground state of H c is a chiral symmetry breaking Wigner lattice 
where either the even or odd sublattice is completely occupied and the other sublattice 
is completely empty. Thus, we expect chiral symmetry breaking in this case too, with 
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critical behavior in the universality class of the Ising model. However, we point out that 
the vacuum energy is not a smooth function of N, being of order e 2 when N is odd and of 
order 1/e 2 when N is even. 

In 2+1-dimensions, our results suggest chiral symmetry breaking, at least for N=l 
and 2, in agreement with the analysis of quenched Schwinger-Dyson equations [16]. It is 
also known that at large N other phases, such as the flux phase itself [17] compete with 
the Neel ordered phase of the antiferromagnet. It would be interesting to examine this 
further, in particular to obtain an estimate of upper critical N, if one exists, in the context 
of the present model. 

Salmhofer and Seiler [8] have given a proof that QED defined on a spacetime lattice 
with > 4 flavors and in > 4 dimensions breaks chiral symmetry. Their model differs 
from ours in the order of symmetry of the strong coupling versus the weak coupling limit. 
Their 4-flavor QED has only Z 2 symmetry at strong coupling and no Goldstone bosons. 
Furthermore, the order parameters differ in the two cases, theirs being in the chiral U(l) 
channel. It is also a puzzle to us that the vacuum energy of their model in the large e 2 
limit is of order one, whereas we find that it is always either of order e 2 or 1/ e 2 depending 
on whether N is even or odd. Resolution of these differences is an important problem. 
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